Efficiency of heat engines coupled to nonequilibrium reservoirs 
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We consider quantum heat engines that operate between nonequilibrium stationary reservoirs. 
We evaluate their maximum efficiency allowed by the second law of thermodynamics. We show that 
it can be expressed in terms of an effective temperature that depends on the nature of the reservoirs. 
We further compute the efficiency at maximum power for different kinds of engineered reservoirs. 
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Engines are devices that convert various forms of en- 
ergy into useful mechanical work and motion. In ther- 
modynamics, two different kinds of machines can be dis- 
tinguished. On the one hand, there are heat engines 
that operate between two reservoirs at different temper- 
ature, such as internal combustion engines [TJ[2]. On the 
other hand, there are molecular motors that are driven 
from equilibrium by varying external parameters, while 
in contact with a single isothermal reservoir [3] |4j . The 
latter describe biological motor proteins as well as arti- 
ficial nanomachines [5] . An essential characteristic of 
any machine is its efficiency defined as the ratio of work 
output to energy input. Whereas for heat engines the ef- 
ficiency is limited by the Carnot formula, rj c = 1 — T1/T2, 
where T\ and T2 are the temperatures of the two ther- 
mal reservoirs (T\ < T%) , it can reach unity for molecular 
motors [7H5]. Maximum efficiency usually corresponds 
to quasistatic conditions, and therefore to zero power. A 
practically more relevant quantity is thus the efficiency at 
maximum power which for heat engines is given by rj c /2 
for small temperature differences [T0TU2| . For molecular 
motors, the efficiency at maximum power can reach the 
thermodynamic limit 1 for strong driving [131 114) . 

Heat engines are usually assumed to be in contact with 
two equilibrium reservoirs. In this paper, we investigate 
the more general case where the engine runs between 
stationary nonequilibrium reservoirs. In a sense, this sit- 
uation interpolates between traditional heat engines and 
molecular motors. Indeed, the efficiency of these heat en- 
gines can be larger than the Carnot efficiency and they 
can operate isothermally. Our study is motivated by the 
recent advent of reservoir engineering techniques in quan- 
tum optical systems, such as ion traps P~5][T|)], microwave 
cavities [XT] [TS] , optical lattices [TH] [20] and optomechan- 
ical systems [5T] , that enable the preparation of nonther- 
mal environments. In addition, theoretical studies have 
shown in individual cases that the efficiency of heat en- 
gines coupled to nonthermal quantum coherent [22] or 
quantum correlated [23] reservoirs can sometimes exceed 
the Carnot limit. The two fundamental questions that 
we here address are therefore: i) under what conditions 
can the efficiency be enhanced beyond Carnot, and ii) 
what is the maximum (universal) efficiency that can be 
reached? In the following, we consider a quantum heat 
engine coupled to general stationary nonthermal reser- 



voirs. It will be convenient to regard these reservoirs 
as perturbed thermal reservoirs, although the perturba- 
tion need not be small. They will then be characterized 
by a temperature and a second parameter (or more) that 
quantifies the deviation from equilibrium, such as the de- 
gree of quantum coherence |22j or the amount of quan- 
tum correlations [53] . We begin by performing a detailed 
analysis of the quantum Otto cycle for a time-dependent 
harmonic oscillator, a paradigm of quantum heat engines 
and a generalization of the common four-stroke car en- 
gine [2~4TI2"5] . We evaluate its efficiency which we express 
in terms of the Hamiltonian of mean force [29] [30], a 
quantum extension of the potential of mean force well- 
known in the statistical theory of fluids [51] . We further 
derive an explicit expression for the maximum efficiency 
of any heat engine allowed by the second law of thermo- 
dynamics; this efficiency may be smaller or larger than 
the Carnot efficiency depending on the properties of the 
reservoirs. We will use the latter to obtain a criterion for 
the increase of the efficiency beyond the Carnot bound. 
Finally, we compute the efficiency at maximum power 
of the Otto engine for the concrete examples of strongly 
coupled, coherent and correlated quantum reservoirs. 

Quantum Otto engine. We consider a quantum Otto 
engine whose working medium is a harmonic oscillator 
with time-dependent frequency oj t [24- 2SJ . The Otto cy- 
cle consists of two isentropic processes during which the 
frequency is unitarily varied between u>i and U2, and of 
two isobaric (constant frequency) processes during which 
the oscillator is connected to two different reservoirs (see 
Fig. 1). A concrete scheme to experimentally realize such 
an engine using a single ion in a linear Paul trap has been 
proposed in Ref. [55] • In the usual Otto cycle, the two 
reservoirs are assumed to be thermal and characterized 
by the inverse temperatures /3.; = l/^^X^), (i = 1,2), 
where ks is the Boltzmann constant. Here, we examine 
the situation where the engine is alternatingly coupled 
to engineered nonthermal reservoirs. For simplicity, we 
will first focus on the case where only the hot reservoir is 
nonthermal. When connected to this reservoir, the oscil- 
lator relaxes to a nonequilibrium state that we write in 
the general form, P2 = cxp[— ^2(^2 + Af/2)]/ '^i where 
H2 = p 2 /(2m) + mw|x 2 /2 is the Hamiltonian of the os- 
cillator at frequency u>2 (m denotes the mass) and 
the normalization constant. The operator A.f/2 quanti- 
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FIG. 1: (Color online) Energy-frequency diagram of the quan- 
tum Otto engine. The thermodynamic cycle consists of two 
isentropic (stroke 1 and 3) and two isochoric (stroke 2 and 
3) processes. In the latter, the heat engine is coupled to en- 
gineered nonthermal stationary reservoirs that are described 
by a temperature and additional parameters that characterize 
the deviation from thermal equilibrium. 



fies the departure from equilibrium and may be arbitrary 
[32] . We note that the operator H 2 = H 2 + A7J 2 can be 
seen as a Hamiltonian of mean force [2^1 [3D] . It is im- 
portant to stress that the state p 2 should be centered, 
(x) = (p) = 0, so that the nonequilibrium reservoir is 
a proper heat source that can only exchange heat with 
the engine, but no work [53]. The mean energy of the 
harmonic oscillator in the nonequilibrium state is then, 



(H*) = Huj 2 (fi 2 + An +1/2), 



(1) 



where ni = [exp(h/3iUJi) — is the mean occupation 
number of a thermal quantum oscillator and An that 
associated with the deviation from the thermal state. 

During the two isentropic parts of the thermodynamic 
cycle (stroke 1 and 3), the time-dependent oscillator is 
isolated and its dynamics is hence unitary. As a result, its 
Schrodinger equation can be solved exactly using a Gaus- 
sian wave function ansatz and the mean energy can be 
evaluated analytically [34 36 . During the two isochoric 
branches (stroke 2 and 4), the oscillator relaxes respec- 
tively to a nonequilibrium and an equilibrium state. The 
corresponding average energies of the harmonic oscillator 
at the four corners of the quantum Otto cycle are, 



(H) 
(H) 
(H) 

(H) n = 



= tvjJi{n x + 1/2), (2a) 

= foj2QJ(ni + l/2), (2b) 

= huo 2 {n 2 + An +1/2), (2c) 

= hw 1 Q*(n 2 + An + 1/2), (2d) 



where we have used Eq. ^ and Eq. (5.12) of Ref. [34]. 
We note that the energies at point C and D are modified 
by the presence of the nonthermal reservoir. The two 



parameters Q\ and Q 2 characterize the degree of adia- 
baticity of the compression and expansion phases 1 and 
3 [34] . Their explicit expressions for any given frequency 
modulation u t , can be found in Refs. [351 136| . They are, 
for example, equal to one for adiabatic processes and to 
(uif + oj 2 ) / (2luiuj 2 ) for a sudden frequency change. 

To evaluate the efficiency of the Otto engine, we need 
to compute work and heat along the four branches of the 
cycle (see Fig. 1). The mean works, denoted by (Wi) and 
(W 3 ) , done during stroke 1 and 3 are given by, 

(Wi) - (H) B -(H) A = (hu 2 Q* 1 -hL> 1 )(n 1 + ±), (3) 
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At the same time, the mean heats, (Q 2 ) and (Qf), ex- 
changed with the reservoirs during stroke 2 and 4 read, 



(Qa) 



= (H) C -(H) B 

= hw 2 (n 2 +An+^)-hoj 2 Q* 1 {n 1 + ^), (4) 
= (H) A -(H) D 

= frwi(fii + -) - hujiQ 2 (n 2 + An + -). 



The efficiency, defined as the ratio of the total work per 
cycle, (W) = — ((Wi) + (W3)), to the heat received from 
the hot reservoir, (Q 2 ), can then be written as, 



n 



Uh (ni + 1/2) - Q* 2 {n 2 + An + 1/2) 
w 2 Ql (fii + 1/2) - (n 2 + An + 1/2) ' 



(5) 



The above quantum expression is exact. It gives the finite 
time efficiency of the quantum Otto engine for any fre- 
quency modulation cj f , any inverse temperature Pi, and 
any nonequilibrium stationary reservoir. 

Maximum efficiency. The second law of thermody- 
namics limits the efficiency of any heat engine, in contact 
with either equilibrium or nonequilibrium reservoirs. We 
next derive the maximum efficiency of a heat engine oper- 
ating between a cold thermal reservoir and a hot nonther- 
mal reservoir using Eq. |5]). According to the Clausius 
statement of the second law "Heat cannot, by itself, pass 
from a colder to a hotter body" [37] . As a result, heat 
is absorbed from the hot reservoir, (Q 2 ) > 0, and flows 
into the cold reservoir, (Qf) < (33J. The latter yields 
the following two conditions: 



n 2 + An+ 1/2 
ni + 1/2 



> 



ni + 1/2 
n 2 + An + 1/2 



< 



(G) 



Maximum efficiency is achieved for adiabatic frequency 
modulations which correspond to Q* 2 = 1. Combining 
Eqs. ([5| and ^ in the high-temperature limit, n^ + 1/2 ~ 
l/(fi,/3jCLij), we therefore find, 

r 1 <n c + h(3 2 uj 1 Ah = n c + p 2 —(AH) c , (7) 

u> 2 
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where (AH) C = (H) c - (H) c £ = huj 2 (An) is the devi- 
ation of the energy of the oscillator from its equilibrium 
value at point C. Equation Q reduces to the Carnot 
efficiency, the universal upper bound for any heat engine 
running between thermal reservoirs when (AH) C = 0. In 
the general case, it extends the Carnot formula to heat 
engines operating between a thermal and a nonthermal 
stationary reservoir. Efficiency nf\ exceeds the Carnot 
efficiency r\ c if and only if (AH) C > 0. The latter condi- 
tion corresponds to a larger area of the thermodynamic 
cycle Fig. 1, and thus to a larger work output. 

Expression Q agrees with the maximum efficiency of 
a quantum Carnot engine in contact with either a quan- 
tum coherent [22j or a quantum correlated |23j reservoir 
(see below) , obtained in the limit of small deviations from 
equilibrium. Equation Q, however, is not limited to the 
close to equilibrium regime; it thus provides an extension 
of the results derived in Refs. [22j [23] to far from equilib- 
rium reservoirs. Moreover, using Eq. (J6|, we may rewrite 
the maximum efficiency |7| in the form, 



= 1 - 



/3i(l + fa (AH} C ) ' 



(8) 



The latter only depends on the two temperatures and on 
the energy deviation from equilibrium at point C. We 
observe that the maximum efficiency is smaller than 1 
and that the nonthermal reservoir can be described by 
an effective temperature, T% s = T% + (AH) C /fcs, that 
may be larger or smaller than that of the unperturbed 
thermal reservoir. This situation is reminiscent of that of 
molecular motors where the external driving is sometimes 
regarded as a nonequilibrium reservoir with an effective 
temperature (see e.g. Ref. [4], Sect. 3.4.2). We note, 
in addition, that the fact that Eq. ^ appears as the 
maximum efficiency of different kinds of quantum heat 
engines (Otto and in particular Carnot) strongly hints at 
its universal validity for all heat engines. 

An important observation is that in thermodynam- 
ics the two thermal reservoirs are supposed to be given 
PJ[2]. In particular, the energetic cost of preparing, say, 
a high temperature reservoir in addition to an ambient 
low temperature reservoir is not taken into account in 
the calculation of the efficiency of a heat engine |55] . 
Such an inclusion would indeed lead to vanishing effi- 
ciencies due to the large (strictly speaking infinite) en- 
ergy content of a proper heat reservoir. We here follow 
the same approach and consider the nonthermal reser- 
voir as given. In this framework, the Carnot formula 
appears as an expression of the second law of thermo- 
dynamics for a particular form of nonequilibrium (two 
thermal reservoirs at different temperature), whereas the 
efficiency ([8| applies to a more general form of nonequi- 
librium (one thermal and one nonthermal reservoir). 
In a similar manner, the maximum efficiency of 1 for 
molecular motors is a consequence of the second law 
for this yet different type of nonequilibrium |39j . It 
is worth noticing that Eq. Q yields a non-zero result, 
?7max = /?2 (AH) C /(l + /3 2 (AH) C ), for isothermal reser- 



voirs, Pi = /?2- In this situation, which is akin to that of 
molecular motors, the nonthermal reservoir can be seen 
as an external nonequilibrium driving. 

Efficiency at maximum power. The efficiency at maxi- 
mum power is often a more relevant quantity than the 
maximum efficiency which corresponds to zero power 
|40) . In contrast to the latter, however, there does not 
seem to be a universal expression for the efficiency at 
maximum power; it not only usually depends on the op- 
timization procedure, but also, as we will show, on the 
details of the nonthermal reservoir. The power output 
of a heat engine is defined as P — (W) /r, where r is 
the duration of the cycle. In the following, we evaluate 
the maximum efficiency and the efficiency at maximum 
power for three different examples. For simplicity, we fo- 
cus on adiabatic compression and expansion, Q\ 2 = 1, 
since nonadiabatic processes lead to smaller efficiencies 
[26l[28], and to the high-temperature regime fahwi <C 1. 

(1) Let us begin with a heat engine strongly coupled 
to a hot reservoir. Strong coupling has recently been 
achieved experimentally in circuit QED |41j and optical 
systems [42] . In standard thermodynamics, the coupling 
to the reservoir is assumed to be vanishingly small |43j . 
For strong coupling, the stationary state of the quantum 
oscillator departs from its equilibrium expression and is 
characterized by the temperature of the reservoir and by 
the coupling strength 7. The corresponding deviation of 
the mean occupation number is given by |30j . 



An, 



24 oj 2 



(9) 



where loo is the (large) cutoff frequency of the reservoir, 
modeled as an ensemble of harmonic oscillators. Since 
An 7 > 0, the maximum efficiency Q of the heat en- 
gine is larger than the Carnot efficiency. The total work 
produced by the engine is, 



-(W) = 



- 1 



1 

W2 



H 2 fajLJD 



UJl 
U>2 



- 1 



(10) 

Assuming that the initial frequency of the oscillator u>i 
(as well as 7, fa ,/3 2 and the cycle time) are fixed and 
by optimizing with respect to the second frequency aj 2 , 
we find that the power is maximum when lui/uj 2 = 

+ ^ 2 /^27 w -d/ 24 ])- As a result, the efficiency 
at maximum power is given by, 



''h 



fa (I + fa (AH) c ) 



(11) 



Equation (11) reduces to the Curzon-Ahlborn expres- 
sion [40) for vanishing coupling and generally exceeds 
it when (AH) C > 0, that is, 7 > 0. Remarkably, 
Eq. (11) can be expressed in terms of the same effec- 
tive temperature as Eq. In the limit of small tem- 
perature differences and small (AH) C , we moreover have 
?7 7 ~ ?? c /2+/3| (AH) C /(2fa) > rj c /2; this result therefore 
lies beyond the range of the usual linear regime [Hfl - fT^] . 
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(2) We now consider a quantum photo-engine made 
of a single mode (the harmonic oscillator) in a resonant 
cavity with moving mirrors, and coupled to a beam of 
thermal two-level atoms that pass through the cavity 
[2"2l 125] . When the atoms are uncorrelated, the beam 
plays the role of a thermal reservoir. By contrast, for 
correlated atoms the engineered reservoir is nonthermal. 
For pairwise thermally entangled atoms, the deviation of 
the mean occupation number from equilibrium is [23] , 



3 2 ftA 2 
4w 2 



AnA:2 



A 



(12) 



when respectively one or the two atoms of a correlated 
pair fly through the cavity, in the limit of high tempera- 
ture, PihiOi <C 1, and weak correlation, ftTiA <C 1. Here 
A is the strength of the interaction that created the ther- 
mal entangled pair. In both cases, the deviation of the 
mean occupation number is inversely proportional to the 
frequency w 2 , as for the strongly coupled reservoir dis- 
cussed above. As a consequence, the efficiency at maxi- 
mum power of the heat engine is given by Eq. ( 11 ). We 
note that the maximum efficiency ^ is larger than the 
Carnot expression when one atom flies through the cavity 
and smaller when the two atoms of a pair pass through it. 
Likewise, the efficiency at maximum power (11) is larger 



or smaller than the Curzon-Ahlborn efficiency whether 
one or two atoms fly through the cavity. 

(3) A quantum coherent nonthermal reservoir can be 
created by sending through the optical cavity a beam 
of thermal three-level atoms whose degenerate ground 
states are prepared in a coherent superposition with rel- 
ative phase (j>. In the high-temperature limit, the devia- 
tion of the mean occupation number is |22) , 



An^ 



1 



(/t/3 2 w 2 ) : 



e cos i 



(13) 



where e is proportional to the amplitude of the atomic co- 
herence. Here An is inversely proportional to the square 
of the frequency uj 2 . The total work produced by the 
heat engine can be readily written as, 



(W) = 



1 



UJ 2 



Pi w 



- 1 



1 



UJ 2 



-11- 



ecoscp 

hf3 2 UJ 2 



(14) 

By maximizing the power with respect to w 2 , keeping 
all other parameters constant as before, we find that the 
power is maximum when, 




1 - 



ecoscp 



1 - 2 




(15) 



in a perturbation expansion for small values of e 
The resulting efficiency at maximum power is then 



V<j> = i - 



ecoscp 
2h(3 2 uj 1 




(16) 



The maximum efficiency Q and the efficiency at maxi- 
mum power ( 16 1 exceed their thermal counterparts when 
the condition cos cj) < is satisfied. 

Generalization. The above results can be extended to 
situations where the two reservoirs are nonthermal. For 
instance, the efficiency of the quantum Otto engine is, 



77 = 1- 



wi (rli + Ani + 1/2) - Q* 2 {n 2 + An 2 + 1/2) 



uj 2 (rli + Ani + 1/2)Q* - (n 2 + An 2 + 1/2) ' 

(17) 

where Ani is the deviation of the mean occupation num- 
ber for the nonthermal reservoir i (i = 1,2). The corre- 
sponds high-tcmperature maximum efficiency follows as, 



W = Vc + ( — (Afl) c - (AH) A ) , (18) 
,w 2 



fa 



(3 1 [l + p 2 ((AH) c -(u J2 /io 1 )(AH) A Y 

where (AH) A = (H) A — (H) A q is the deviation of the en- 
ergy of the oscillator from its equilibrium value at point 
A. The efficiency ( fl8| ) is larger than the Carnot effi- 
ciency when lji (AHjg — lo 2 (AH) a > 0. It is inter- 
esting to note that Eq. (18) can again be expressed in 
terms of an effective temperature, T| ff = T 2 + ((AH) C — 
(cj 2 /oj 1 ) (AH) A )/kB, that only depends on the average 
deviations at point A and C and not at points B or D. 

So far we have implicitly assumed that the reservoir 
with the highest temperature also leads to the oscillator 
state with the highest energy, and vice versa. While this 
is true for thermal reservoirs, it is not necessarily true 
for general nonthermal reservoirs; the average energy at 
point A could, for example, be higher than the average 
energy at point C . This issue can be addressed by adapt- 
ing the definition of the efficiency ^ and generalizing the 
Clausius statement of the second law, replacing hot by 
high energy and cold by low energy. 

Conclusions. Combining the efficiency of a quantum 
Otto engine and the second of law of thermodynamics, 
we have computed the maximum efficiency of a heat en- 
gine operating between nonthermal stationary reservoirs. 
We have shown that the latter can be expressed in terms 
of an effective temperature which characterizes the devi- 
ation from thermal equilibrium. We have further estab- 
lished that the maximum efficiency exceeds the Carnot 
bound when the effective temperature is larger than the 
high temperature. We have in addition evaluated the ef- 
ficiency at maximum power for three cases of interest. 
Our results provide a theoretical framework for a new 
class of engineered heat engines that interpolate between 
standard heat engines and molecular motors. 
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